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ABSTRACT 

A meta-analysis is performed of the literature on evolution in cosmic star-formation 
rate density from redshift unity to the present day. The measurements are extremely 
diverse, including radio, infrared, and ultraviolet broad-band photometric indicators, 
and visible and near- ultraviolet line-emission indicators. Although there is large scatter 
among indicators at any given redshift, virtually all studies find a significant decrease 
from redshift unity to the present day. This is the most heterogeneously confirmed 
result in the study of galaxy evolution. When comoving star-formation rate density is 
treated as being proportional to (1 + z)^ , the meta-analysis gives a best-fit exponent 
and conservative confidence interval of /? = 2.7 it 0.7 in a world model with (0^/ , Oa) = 
(0.3,0.7) and /3 = 3.3 ± 0.8 in {Qm,^a) = (1.0,0.0). In either case these evolutionary 
trends are strong enough that the bulk of the stellar mass at the present day ought to 
be in old (> 6 Gyr) populations. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: stel- 
lar content — history and philosophy of astronomy — methods: statistical — stars: 
formation 

1. Introduction 

The study of galaxy evolution is filled with negative results. Despite intensive efforts, only 
small or subtle evolution has been detected in the number density of field galaxies (eg, Lilly et al. 
1995; Heyl et al. 1997; Hogg 1998; Lin et al. 1999; Cohen 2002), in their masses (eg, Vogt et al. 
1996, 1997; Treu et al. 1999; Brinchmann & Ellis 2000; Cohen 2002), or in their clustering relative 
to "stable clustering" (eg, Le Fevre et al. 1996; Carlberg et al. 1997; Small et al. 1999; Hogg et al. 
2000; Carlberg et al. 2000), from redshift unity to the present day. At the same time, the constraints 
on galaxy evolution have not been made strong enough to allow definitive results from the classical 
cosmological tests. Has the galaxy evolution community got anything positive to say? Indeed it 
has. 



- 2 - 



For a long time it has been observed that apparently faint and (at one time presumed to be, now 
largely known to be) distant galaxies are, on average, bluer in color than their local counterparts 
(eg, Koo & Kron 1992, and references therein). This was attributed to higher star- formation rates 
or younger ages at earlier times. As the data on distant galaxies has improved, this conclusion has 
been bolstered, with photometric and spectroscopic observations spanning the full electromagnetic 
spectrum. 

Often, the literature on galaxy evolution focuses on disagreements between star-formation rate 
measurements, as it should, since such disagreements ought to point to important issues in selection 
effects, dust extinction, the initial mass function of stars, and stellar population syntheses (eg, Cowie 
et al. 1999; Bell & Kennicutt 2001; Hopkins et al. 2001; Sullivan et al. 2001). What is much more 
remarkable than the disagreements between the measurements, however, is the very important 
respect in which virtually all studies of galaxy star-formation rates agree: The star-formation rate 
density in normal galaxies has been declining from redshift unity towards the present day. 

It is of great importance that the star-formation measurements span a number of different ob- 
servational techniques, a number of different regions in the electromagnetic spectrum, and a number 
of physically different star-formation indicators. The different studies suffer from different selection 
effects and have different sensitivities to dust extinction and the stellar initial mass function. For 
example, as the absorption due to dust in a star-forming galaxy increases, the ultraviolet luminosity 
decreases, the optical Ha line decreases by less, the radio emission from supernovae is unaffected, 
and the far-infrared emission actually increases. The different measurements are also performed 
by different groups of investigators, with different approaches and different preconceptions. This 
heterogeneity of the confirming evidence makes the decline in star-formation rate density the most 
secure result in galaxy evolution, and one of the most secure in all of astrophysics. 

In what follows, a quantitative meta-analysis of the star-formation history literature is per- 
formed, in which the declines inferred from different observational studies are compared in a fair 
way, to establish the consistency of the literature and combine the measurements responsibly. Be- 
cause of the "responsibility" constraint, many studies relevant to questions of star- formation rate 
have been dropped from the meta-analysis; even though almost all of them support the final con- 
clusions. 

The star-formation rate measurements in question are measurements of mean comoving cosmic 
star formation rate density p*, the average mass in stars formed per unit time per unit comoving vol- 
ume. As discussed below, except where specifically noted, all reported measurements are corrected 
to a cosmological world model with (/i, Om,^^a) = (0.7,0.3,0.7), where Hq = 100 h kms^^Mpc"^ 
and Qm and are the present-day, scaled densities of matter and cosmological constant (eg, Hogg 
1999). 
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2. 



Star formation indicators 



As emphasized above, confidence in the result that star formation is evolving comes from from 
the heterogeneity of star formation indicators that show similar evolution. 



Short-lived stars tend to have high temperatures, so the ultraviolet continuum of galaxies is an 
indicator of recent star formation. With optical telescopes, near-ultraviolet wavelengths A ~ 2500 A 

can be observed in the U (or similar) bandpass at redshifts z > 0.4. A stellar population created 
in a short burst is bright at 2500 A for 20 Myr (eg, Bruzual A. & Chariot 1993; Leitherer & 
Heckman 1995) because only short-lived stars emit at these wavelengths. 

The near-ultraviolct luminosity of a galaxy is a direct measure of its young stellar population, 
and therefore a good indicator of recent star formation activity. The primary disadvantage of this 
indicator is that it is strongly affected by dust extinction; correction for extinction requires good 
knowledge of the reddening law in the near-ultraviolet and the geometry of the dust relative to the 
stars (eg, Calzetti et al. 1994). Another disadvantage is that ultraviolet continuum can be produced 
by nuclear activity; in fact possible contamination by nuclear activity is a disadvantage of all of 
the indicators. 

At redshifts z < 0.4, the 2500 A emission can only be estimated from the ground by extrap- 
olation from the visual and near- ultraviolet. This is a disadvantage, because the extrapolation 
depends on the shape of the galaxy spectral energy distribution in the ultraviolet; this is very 
sensitive to small amounts of dust or small changes in star formation history and activity. The 
forthcoming GALEX mission (Martin et al. 1997), which will obtain ultraviolet imaging of normal 
nearby galaxies, will solve this problem. 



The Ha line is emitted in the recombination spectrum of hydrogen; it is a star formation 
indicator because only short-lived stars produce significant luminosity in photons with energies 
above the hydrogen ionization energy 13.6 eV. 

Because it involves ionization and recombination, the Ha line is in some sense a less-direct 

indicator of recent star formation. However, it is weighted towards even higher-mass (and therefore 
shorter-lived) stars. It is also less affected by dust extinction (though not completely insensitive). 
The indicator has the disadvantage that in a significant fraction of galaxies, the total Ha luminosity 
can be dominated by an active nucleus. 

At redshifts near unity. Ha can only be measured with near-infrared spectrographs. Some such 



2.1. 



Rest-frame ultraviolet continuum 



2.2. 



Ha emission line 
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instruments are now on-line, but measurements of line strengths in faint galaxies have not been 
made in large numbers. This is likely to change in the near future. 

2.3. Near-ultraviolet emission lines 

Ultraviolet radiation from stars ionizes many atomic species, not just hydrogen, and several 
of these species emit bright lines in the optical and near-ultraviolet in their recombination spectra. 
The most prominent and useful line is the forbidden [O II] line at 3727 A, because it is bright and 
visible in optical spectrographs in a wide redshift range 0.2 < z < 1.3. Because the line is forbidden, 
it is a collisionally excited transition emitted by Oxygen in ionization regions with electron densities 
n ^ 10^ cm~^. These conditions are common around very young and forming stars. 

In the local Universe, [O II] emission appears to be closely related to Ha emission (Kennicutt 
1992), so it is a reasonable tracer of star formation. However, its strength depends on the metallicity, 
density, and dust content of the nebular gas, so the relationship between line luminosity and star 
formation rate is expected to have a much greater variance than that for the Ha line. Its primary 

advantage is in redshift coverage. 

The [O II] line shares the disadvantage that it can be contaminated by nuclear activity. 

2.4. Far-infrared continuum 

Star forming regions are often very dusty; thick dust surrounding young stars will absorb 
ultraviolet and optical emission and re-emit it at far-infrared wavelengths 30 < A < 200 /xm. Since 
short-lived stars are far more luminous than long-lived stars (by a factor which overwhelms their 
relative scarcity in typical initial mass functions), the far- infrared luminosity of a galaxy may be 
approximately the bolometric luminosity of its dust-enshrouded, short-lived stars. This makes it a 
star-formation indicator. 

Far-infrared emission is a less direct indicator than ultraviolet flux, but it has opposite sen- 
sitivity to dust; dustier galaxies are better measured, while pristine galaxies arc missed. It has 
the enormous disadvantage that the far infrared cannot be measured directly at ground-based tele- 
scopes. It must be inferred from sub-millimeter and near-infrared measurements. The study of 
far-infrared emission will be revolutionized with the SIRTF mission. There is an additional prob- 
lem that the far-infrared emission can have a substantial contribution from an active nucleus, and 
the visual extinction makes it difficult to differentiate sources of luminosity (eg Sanders &; Mirabel 
1996, and references therein). 
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2.5. Radio continuum 

A short time ~ 30 Myr after a burst of star formation, stars with masses M > 8 Mq complete 
most of their nuclear burning and explode in what are expected to be Type II and Type lb super- 
novae, although the evidence for this association is not iron-clad (eg, Filippenko 1997, and references 
therein). These supernovae and their remnants are bright at radio frequencies 1 < v < 100 GHz 
for t ~ 100 Myr (eg, Condon 1992, and references therein). 

Radio observations have the great advantage that they axe insensitive to dust. Star formation 

measurement is not rocket science, but the radio indicator is the least well-calibrated, simply 
because the energetics of radio supernovae are less well-understood than other aspects of stellar 
populations. Like the other indicators, radio measurements can be contaminated by active nuclei, 
which, at faint levels, are not easily distinguished from supernovae. 

2.6. Other methods 

Many other methods have been suggested but wait for new data. It has been suggested that 
the X-ray sources created by stellar evolution could be used as an indicator (Ghosh k, White 2001). 
Much work is going into understanding the winds from forming and young stellar populations (eg, 
Kudritzki & Puis 2000. and references therein); one can imagine measuring star formation rates 
using mechanical, rather than electromagnetic, luminosity. When the SIRTF mission is able to 
calibrate the relationships between strong mid-infrared spectral features from interstellar molecules 
and the ultraviolet radiation fields that excite the transitions (eg, Li & Draine 2002) , there will be 
a new industry of understanding star formation activity via the mid-infrared. 

3. Method and results 

3.1. Studies excluded 

This meta-analysis is of star-formation rate variation within a single observational study, at 
redshifts z < 1. For this reason, no study was included in this meta-analysis if it did not report 
more than one independent z < 1 star-formation rate measurement. This excluded several otherwise 
relevant studies (eg, Connolly et al. 1997; Madau et al. 1998; Treyer et al. 1998; Yan et al. 1999; 
Sullivan et al. 2000; Thompson et al. 2001). An exception was made for a group of very similar 
surveys for Ha line luminosity density, which were put together to make a "combined Ha" study, 
described below. 

Several relevant studies were dropped because the targets were selected or the results were 
presented such that the results could not be fairly used as cosmic star-formation rate measurements 
(eg, Schade et al. 1996; Cowie et al. 1997; Guzman et al. 1997; Lilly et al. 1998; Blain et al. 1999), 
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or because estimates of measurement uncertainties were not provided (eg, Cram 1998). 

Only publications appearing in the refereed literature prior to 2002 August 1 were considered. 



The cuts left ten studies (Lilly et al. 1996; Hammer et al. 1997; Rowan-Robinson et al. 1997; 
Hogg et al. 1998; Cowie et al. 1999; Flores et al. 1999; Mobasher et al. 1999; Haarsma et al. 2000; 
Jones &; Bland-Hawthorn 2001). Notes on individual studies are given in Table 1, but in most 
cases, the confidence interval on each star-formation-ratc-dcnsity p^, point was measured (with a 
ruler) from the relevant Figure in each paper. The measured points are shown in Figure 1. At any 
redshift, there is a large scatter. 

In addition to the primary ten studies, three very similar studies of the cosmic Ha line lumi- 
nosity density (Gallego ct al. 1995; Tresse & Maddox 1998; Glazebrook ct al. 1999), each measuring 
the density at a different redshift, were combined into a single "combined Ha" study, which was 
included as if it were an additional individual study. 



Each data point, centered on some redshift z, was crudely "corrected" to the default world 
model used in this paper by the correction factor 



where Dl and Vc are the luminosity distance and comoving volume (eg Hogg 1999) out to z in the 
cosmological world model employed in the study in question, and Dl and Vc are the same but for 
the fiducial world model (/i, VLm, ^k) = (0.7, 0.3, 0.7) used here. This world-model correction makes 
two assumptions: The first is that the star-formation indicators have been calibrated somehow 
independently of Hubble constant, in the sense that a change in the Hubble constant changes the 
inferred luminosity and inferred star formation rate by the same factor. The second is that each 
reported data point can be treated as well-localized at its redshift; for the precisions reported in 
these studies this ought to be fine. 

Linear fitting to the model p* oc (1 -|- z)l^ was performed in log(/9*) space with points weighted 
according to their reported uncertainties. By this procedure, the confidence intervals reported in 
the individual studies are treated as representing Gaussian uncertainties in log(/9*) space. In detail, 
of course, this assumption is incorrect, but since none of the studies is extremely precise and since 
most of the reported uncertainties are not accurate (most only include the Poisson contribution) , 
the assumption does not significantly affect the results. The best-fit (3 values and the linear-fitting 
uncertainties provided by the covariance matrix are given in Table 1 and shown in Figure 2. 



3.2. 



Studies included 



3.3. Data fitting 




(1) 
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3.4. Results 

The weighted mean (using the inverse variances from the hnear fits as weights) of the best-fit 
P values is /? = 2.74 it 0.28, a nominal ten-sigma result. 

All but two of the /3 measurements are consistent (at ~ 2 a) with the weighted mean. In detail, 
the weighted mean relies on the assumption that the studies have reported accurate uncertainties; 
but weighted and unweighted mean values barely differ. Of the two outliers, one (Hammer et al. 
1997) appears to be an outlier simply because it carries such a small uncertainty; most of the 
studies report only Poisson contributions to their error budgets, so the uncertainties are better 
treated as lower limits. The other outlier (Mobasher et al. 1999) is more problematic, but the 
sample is small (only 23 galaxies in the z ~ 0.5 point), the sample is cut on both radio power and 
optical magnitude (at a fairly bright R < 21.5 mag), and sources spectrally classified as having 
active nuclei were excluded (roughly 25 percent of the sample). 

Given the underestimated measurement uncertainties, a much more conservative method for 
computing a mean with a reasonable confidence interval is the bootstrap resampling technique. In 
each of 10^ trials, ten studies were randomly chosen (with replacement) from the group of ten and 
a weighted mean (using inverse variances from the linear fits as weights) for (3 was computed. Note 
that each trial will have some studies repeated, some missing. From the 10^ weighted means, the 
lowest and highest 16 percent were discarded, leaving the central 68-percent confidence region. This 
region has a central value and extent of /3 = 2.7 ± 0.7. This is a conservative but robust estimate 
of the average value for p. As described above, all measurements have been corrected to a world 
model with (J^m,J^a) = (0.3,0.7). 

If all results are transformed to an Einstein-de Sitter Universe with {^m,^a) = (1.0,0.0), the 
bootstrap-resampling 68-percent confidence interval is /? = 3.3 it 0.8. 

4. Discussion 

Despite the scatter in Figure 1, the decline in comoving star-formation rate density p* from 
redshift unity to the present day is well confirmed by multiple studies. The best average value 
of the evolution exponent /3, in the parameterization oc (1 -|- z)^ , is /3 = 2.7 it 0.7. This value 
is only four standard deviations from zero, but the confidence interval is computed in the most 
conservative way (resampling); the straightforward weighted mean of the measurements gives a 
ten-sigma result. 

In addition, almost all faint-galaxy studies which are relevant, but for one reason or another 
could not be included in this meta-analysis, bolster the conclusion that the star-formation rate has 
been dropping since redshift unity. Faint galaxies are bluer than nearby galaxies (eg, Koo & Kron 
1992; Smail et al. 1995; Williams et al. 1996). The surface brightnesses, colors, and line emission 
of normal galaxy disks appear to have all been dropping (Schade et al. 1996; Lilly et al. 1998); 
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the same appears true for blue compact galaxies (Guzman et al. 1997). The luminosity functions 
of blue galaxies and star-forming galaxies have been evolving (Lilly et al. 1995; Heyl et al. 1997; 
Cram 1998; Mallen-Ornelas et al. 1999; Lin et al. 1999). Models of mid-infrared and sub-millimeter 
source counts may require a drop in star-formation activity since redshift unity (Blain et al. 1999; 
Elbaz et al. 1999). It is difficult to understand the observed chemical evolution of damped Lya 
absorbers without much higher star formation activity at redshift unity than the present day (Pei 
&; Fall 1995; Pci ct al. 1999). Perhaps more indirect, the fraction of galaxies classified as "irregular" 
appears to be dropping with cosmic time (eg, Griffiths et al. 1994; Abraham et al. 1996; Odewahn 
et al. 1996; van den Bergh et al. 2000), and irregulars, locally, have above-average star- formation 
rates. 

The mere consistency of these star-formation studies is only half the story. It is equally im- 
portant that the studies span a wide range of observational techniques, and a range of physically 
different star-formation indicators. This heterogeneity makes their consistency much more impres- 
sive and much more convincing. In particular the different methods have different sensitivities to 
dust extinction and the slope of the initial mass function of stars; it is hard to argue that the 
observed trends are conspiracies of dust and mass function variability. 

The evolution in star-formation rates is the most secure result in the field of galaxy evolution. 

Incidentally, this also makes it the strongest argument from galaxy evolution against the pure 
steady-state model of cosmic origins (Bondi & Gold 1948; Hoyle 1948). 

If the best-fit value of P and its conservatively estimated uncertainty are taken at face value, 
then P > 1.3 in the (UmMa) = (0.3,0.7) world model and /? > 1.5 in the (J^m,J^a) = (1.0,0.0) 
world model. In cither case, this implies that not only has the comoving star formation rate 
density /9*(-z) been declining since redshift unity, but so has t{z) Pt,{z), the comoving star formation 
density per logarithmic interval in cosmic time. This function t{z) Pt,{z) more accurately represents 
the formation time of the bulk of stars. Because (3 is above the logarithmically divergent value 
(/? > 1.3 in the default world model), most of the stellar mass at the present day ought to be in old 
(> 6 Gyr) stellar populations. This is the primary "prediction" of this meta-analysis; it appears to 
be consistent with the data (Fukugita et al. 1998; Hogg et al. 2002). 

In addition to observations of distant galaxies, we have another powerful and independent 
"fossil record" of galaxy formation and cvohition: Our own Galaxy contains stellar populations 
formed over the entire history of the Universe. The star formation history of the Galaxy can 
be inferred from its internal distribution of apparent stellar ages. Unfortunately, the distribution 
measured in the Solar neighborhood is not consistent with a large decline in the Galaxy's star 
formation rate over the last ~ 8 Gyr (Barry 1988; Pardi & Ferrini 1994; Prantzos & Silk 1998; 
Rocha-Pinto et al. 2000; Gizis et al. 2002). Is the Milky Way atypical? Or, perhaps more likely, 
it may be that measurements of the age distribution in the solar neighborhood not representative 
of the Galaxy as a whole; this is likely if a galaxy's stars are formed from the inside out. The 
resolution of this inconsistency between high-redshift and Milky-Way determinations of the star- 
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formation rate may lead to great progress in cosmology and galaxy evolution. 

It is a pleasure to thank Mike Blanton, Karl Glazebrook, Deborah Haarsma, Heath Jones, 
Jim Peebles, David Schlegel, Rosemary Wyse and Sara Zimmerman for help with the literature, 
scientific discussions, and software. This research made use of the NASA ADS Abstract Service. 
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Table 1. star-formation rate measurements, corrected 



reference 



technique 


z ^ 


logio P* 




Lilly et al. (1996) " 


0.35 


-2.015 ±0.075 


3.93 ± 1.07 


rest-frame ultraviolet 


0.60 


-1.742 ± 0.075 




(0.50,1.00,0.00) 


0.85 


-1.457 ±0.150 




Hammer et al fl997) 


0.48 


-1 797 ± 055 


7.22 ± 1.03 


spectral features, esp. [0 II] 


0.62 


-1.358 ± 0.028 




(0.50,1.00,0.00) 


0.85 


-1.175 ±0.087 




Rowan-Robinson et al. (1997) 


0.55 


-0.815 ± 0.145 


6.93 ± 4.16 


far infrared (inferred) 


0.85 


-0.282 ± 0.285 




(0.50,1.00,0.00) 








Hogg et al. (1998) 


0.20 


-1.827 ± 0.133 


2.88 ± 0.61 


[0 II] line emission 


0.40 


-1.227 ±0.093 




(1.00,0.30,0.00) 


0.60 


-1.141 ±0.067 






0.80 


-1.200 ± 0.078 






1.00 


-0.899 ± 0.097 




Cowie et al. (1999) ^ 


0.37 


-1.890 ±0.122 


1.09 ±1.35 


rest-frame ultraviolet 


0.63 


-1.666 ±0.119 




(0.65,1.00,0.00) 


0.78 


-1.773 ± 0.094 




Flores et al. (1999) 


0.35 


-1.445 ±0.275 


4.47 ± 3.02 


far infrared (inferred) 


0.60 


-1.167 ±0.250 




(0.50,1.00,0.00) 


0.80 


-0.882 ± 0.260 




Mobasher et al. (1999) ^ 


0.15 


-1.553 ±0.080 


-4.05 ± 1.18 


1.4 GHz radio 


0.50 


-2.021 ± 0.110 




(0.50,1.00,0.00) 








Haarsma et al. (2000) 


0.28 


-1.157 ±0.173 


3.61 ± 1.48 


1.4 GHz radio 


0.46 


-0.920 ± 0.171 




(0.50,1.00,0.00) 


0.60 


-0.869 ± 0.209 






0.81 


-0.598 ±0.151 




Jones k Bland-Hawthorn (2001) ^ 


0.08 


-1.892 ±0.301 


2.21 ± 2.76 


Ha line emission 


0.24 


-1.499 ±0.092 




(0.50,1.00,0.00) 


0.40 


-1.501 ±0.166 




Lanzetta et al. (2002) ' 


0.25 


-2.061 ± 0.065 


0.67 ± 0.68 


rest-frame ultraviolet 


0.75 


-1.963 ±0.075 





(1.00,1.00,0.00) 
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Table 1 — Continued 



reference 
technique 




logio P* 


f3' 


combined Ha J 


0.03 


-1.907 ±0.040 


3.75 ±0.51 


Ha line emission 


0.20 


-1.620 ± 0.040 




(0.50,1.00,0.00) 


0.90 


-0.967 ±0.145 




weighted mean 






2.74 ± 0.28 


bootstrap resamphng ^ 






2.68 ± 0.75 



'^Cosmological parameters used in the original reference. For the purposes 
of fitting, the results have been corrected to an effective world model of 
{h,0.M,9\) = (0.7,0.3,0.7); see the text for details. 

is rcdshift. Only points at rcdshifts z < 1 were included in the fits. 

^p^, is the comoving star-formation rate density. All measurements have been 
corrected to the fiducial world model as described in the text. 

is the exponent in parameterization cx (l+z)^. See text for explanation 
of the fitting technique and calculation of uncertainties. 

^Data for Lilly et al. (1996) are taken from analysis of Glazebrook et al. 
(1999). 

^Luminosity density in Cowie et al. (1999) has been converted to star- 
formation rate with the conversion of Cowie et al. (1997). 

^Mobasher et al. (1999) do not give values for the cosmological parameters; 
I have guessed (0.5,1.0,0.0). 

^Integrations of the functions in Jones & Bland- Hawthorn (2001) were kindly 
provided by Jones (private communication). 

'The "scaled break intensity" results of Lanzetta et al. (2002) are used here. 

J The "combined Ha" line is a compilation of three similar studies by Gallego 
et al. (1995), Tresse & Maddox (1998), and Glazebrook et al. (1999), with 
star-formation analysis by Glazebrook et al. (1999). 

'^The central value and uncertainty reported for the bootstrap resampling 
indicate the central 68-percent confidence interval; see the text for details. 
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Fig. 1. — The stax-formation rate measurements used in the meta-analysis, all corrected (crudely; 
see text) to a world model of (/i, Omj^^a) = (0.7,0.3,0.7). See the text and Table 1 for details. 
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Fig. 2. — The best-fit values of evolutionary exponent /? in the parameterization oc (1 + zf for 
the measurements in Table 1 and Figure 1, with world model (P^Mi^a) = (0.3,0.7). The fits are 
performed under a number of assumptions described in the text. The dashed line is the weighted 
mean and the dotted lines indicate the central 68-percent confidence interval from the bootstrap 
resampling described in the text. 



